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Probing d-wave pairing correlations in the pseudogap regime of the cuprate 
superconductors via low-energy states near impurities 
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The issue of probing the pseudogap regime of the cuprate superconductors, specifically with 
regard to the existence and nature of superconducting pairing correlations of d-wave symmetry, is 
explored theoretically. It is shown that if the d-wave correlations believed to describe the super- 
conducting state persist into the pseudogap regime, but with pair-potential phase-fluctuations that 
destroy their long-range nature, then the low-energy quasiparticle states observed near extended 
impurities in the truly superconducting state should also persist as resonances in the pseudogap 
regime. The scattering of quasiparticles by these phase-fluctuations broadens what was (in the 
superconducting state) a sharp peak in the single-particle spectral function at low energy, as we 
demonstrate within the context of a simple model. This peak and its broadening is, in principle, 
accessible via scanning tunneling spectroscopy near extended impurities in the pseudogap regime. 
If so, such experiments would provide a probe of the extent to which d-wave superconducting cor- 
relations persist upon entering the pseudogap regime, thus providing a stringent diagnostic of the 
phase-fluctuation scenario. 
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I. INTRODUCTION 

Among the challenges presented by the high- 
temperature superconductors, one of the most persis- 
tent concerns the anomalous properties of the normal 
state of these materials. In particular, the suppression 
of single-particle spectral weight around the Fermi en- 
ergy [|lril foi' temperatures above the superconducting 
transition temperature of the underdoped cuprates 
indicates that the electronic behavior of these materi- 
als deviates substantially from that of conventional su- 
perconductors. There have been several theoretical sce- 
narios proposed to explain this so-called pseudogap be- 
havior [Q. In the present Paper, we shall be concerned 
with a particular one of these, viz., the phase-fluctuation 
scenario ]lO|-|T3|. According to this scenario, supercon- 
ducting correlations in the form of Cooper pairing are 
presumed to exist (and to be responsible for the loss of 
single-particle spectral weight) for temperatures T below 
the pseudogap onset temperature T* . However, in the 
intermediate temperature range (i.e. < T ^ T*) the 
long-range spatial and temporal coherence in the phase 
of Cooper-pair wave functions (occurring for T < Tc) 
is presumed to present only up to intermediate length 
scales, having been disrupted on longer length scales by 
the presence and motion of vortex excitations. In other 
words, although long-range phase coherence is absent in 
it, the pseudogap regime is quantitatively distinguished 
from the conventional non-superconducting state by the 
presence of substantial, residual, local d-wave pairing cor- 
relations. 

Several recent experimental investigations support the 
notion of the phase-fluctuation scenario as the origin of 
pseudogap phenomenology, including Refs. ||lj,|l^. Fur- 



ther steps towards an understanding the nature of the 
pseudogap regime would be furnished by experimental 
probes that are targeted towards the question of the 
existence of the putative local superconducting correla- 
tions p6|-p^. Such probes would have the potential to 
discriminate between scenarios based on pairing correla- 
tions and those in which the pseudogap is due to some 
other mode of electronic ordering M. The purpose of 
the present Paper is to identify one such probe: scan- 
ning tunneling spectroscopy (STS) measurements of the 
single-particle spectral function near extended impurities 
in the pseudogap regime |20j . 

Before explaining the nature of this probe, let us pause 
to recall one of its essential ingredients. It is has long 
been recognized that the scattering of quasipar- 

ticles between states corresponding to differing signs of 
the d-wave pair-potential leads to the existence of low- 
energy states [|4j . Such scattering, and hence low-energy 
state- formation, occurs, e.g., at suitably oriented sur- 
faces in the cuprates, leading to the observed zero-bias 
anomaly in the tunneling conductance |26|. In the set- 
ting of impurity (rather than boundary) scattering in d- 
wave superconductors, low-energy states, which in this 
case are localized near the impurities, have been ob- 
served in STS experiments |^,^, and have been dis- 
cussed theoretically for the case of point-like and 
extended p9|j3C| ] impurities. We remind the reader that 
these states co-exist with the continuum of low-energy 
quasiparticle states associated with the nodes at which 
the d-wave pair-potential vanishes. Distinguishing be- 
tween impurity states and nodal states is straightforward 
because the former give a peaked contribution to the 
spectral function whereas the contribution from the lat- 
ter vanishes linearly at zero energy. (A clear example of 
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this is furnished by the data reported in Ref. p7[|.) 

Returning to our main task, viz., probing the pseu- 
dogap regime for pairing correlations, we now state the 
central idea on which the present Paper is based. Let us 
suppose that the pseudogap state is indeed distinguished 
by the presence of local (but not long range) d-wave pair- 
ing correlations. The lack of long-range phase coherence 
in such a state complicates the experimental observation 
of the more conventional signatures of superconducting 
correlations (such as the Meissner effect). However, the 
low-energy quasiparticle states occurring near extended 
impurities are localized in space, i.e., their existence only 
relies on the presence of local d-wave pairing correlations. 
In consequence, these states should be only weakly af- 
fected by the destruction of long-range superconducting 
order that occurs at Tc, and hence should "survive" the 
transition into the pseudogap state. Thus, if STS experi- 
ments were to reveal a sharp feature in the single-particle 
spectral function at low voltage bias, this would be ev- 
idence for the presence of local d-wave pairing correla- 
tions. And if such experiments were able to characterize 
the temperature- and doping-dependence of the width of 
this spectral feature, this would provide a characteriza- 
tion of the nature of these finite-range spatial and tempo- 
ral correlations. Of course, the interpretation of any ex- 
periment conducted at nonzero temperature would have 
to contend with broadening (arising, e.g., from thermal 
fluctuations of the sample). In order to minimize the 
consequent smearing of the spectral function, which has 
the potential to obscure the very feature being sought, it 
would be preferable to examine cuprate systems having 
low values of Tc (e.g. heavily underdoped systems). 

Why are we focusing on the case of (non-magnetic) ex- 
tended impurities as a probe of the pseudogap regime? 
The main reason underlying this choice is as follows. In 
the setting of extended impurities (as well as boundaries) 
in the superconducting state, it has been shown that the 
existence of low-energy quasiparticle states is a direct 
consequence of the d-wave nature (or, more precisely, the 
sign-changing nature) of the pair-potential |29) . Extend- 
ing this reasoning to the case of extended impurities in 
the pseudogap regime, there the existence of low-energy 
states strongly depends on the presence of superconduct- 
ing correlations of d-wave type. 

We note that Kruis et al. have studied the den- 
sity of states near a point-like impurity in the context 
of a simple phenomenological picture of the pseudogap 
regime. In this picture, the physics of the pseudogap 
regime is incorporated through the hypothesis that, in 
the absence of the impurity, the single-particle density of 
states vanishes linearly at the Fermi energy. At present, 
the extent to which an approach based on this picture 
can yield information about pairing correlations in the 
pseudogap regime is not clear. 

The present Paper is organized as follows. In Sec. || we 
provide a framework for discussing the influence of local 
d-wave pairing correlations on quasiparticle states near 
extended impurities l33|, focusing on the single-particle 



spectral function. As we shall see, our expression for 
this spectral function will take the form of a density of 
states (determined at a fixed, locally phase-randomized 
d-wave pair-potential) averaged over the fluctuations of 
the phase field. In Sec. Ill we develop a semiclassical 



scheme for computing this density of states at fixed pair- 
potential in which we treat the long-wavelength pair- 
potential phase variations via perturbation theory. This 
scheme allows us to focus on the contribution to the den- 
sity of states at low energies, which involves states arising 
from changes in sign (as the momentum is varied) of the 
local pair-potential. At this point we will have obtained 
an expression for the spectral function which consists of 
a sum of terms each associated with one classical scat- 
tering trajectory that pass es thr ough the vicinity of the 
impurity potential. In Sec. IV A we perform the average 
over the configurations of the fluctuating pair-potential 
arrived at via locally randomizations of the phase, under 
the assumption that the distribution of phase-field config- 
urations is Gaussian. We make an approximation to the 
resulting expression for the spectral function that is valid 
for the case of phase correlations persisting beyond the 
Cooper-pair size, finally arriving at an expression for the 
spectral function near an extended impurity that consists 
of a Gaussian peak the linewidth of which is proportional 
to the typical gradient of the phase of the pair-potential. 
In this section we also compute the linewidth of the spec- 
tral function in the pseudogap regime, doing so by as- 
suming that the phase fluctuations are governed by the 
BKT theory of the two-dimensional ATy-model. finally, 
m Sec. we make a numerical estimate of the linewidth 
by invoking the results of recent high-frequency conduc- 
tivity data 111], and provide some concluding remarks. 



II. MODEL OF CUPRATES WITH AN 
EXTENDED IMPURITY 

In the present section we formulate the task of ob- 
taining the single-particle spectral function in settings 
of systems of fermions interacting via some fermion- 
fermion coupling and also interacting with an exter- 
nal single- fermion potential (which could, e.g., repre- 
sent an impurity potential). As our aim is to address 
the phase-fluctuation picture of the pseudogap regime of 
the cuprates, we envision following the standard field- 
theoretic route (see, e.g., Ref. Q) of exchanging the 
fermion-fermion coupling for a suitable collective quan- 
tum field A. Thus we arrive at the following formula for 
the one-fermion Green function Q{x,x'): 



[ DA^DAe-^^^^ Gix,x';A) 

g{x,x') = L_ 



DA^DAe-'^I^l 



(la) 



dr + h A 
At dr-h 



G{x,x']A) = -5{x-x'), (lb) 
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i.e., a Bogoliubov-de Gennes Green function for pair- 
potential A, suitably averaged over quantum field A, 
the action for which is S[A]. Here, x = (r, r) where 
r and r are respectively the spatial position (in the 
two-dimensional Cu02 plane) and the Matsubara time. 
Furthermore, h = — — kp + V{r), in which fc| is 
the chemical potential [i.e. fcp (= Stt/Af) is the Fermi 
wave vector], V is the single-particle potential, and we 
have adopted units in which /2m = 1, to being the 
(effective) mass of the electrons and holes. The op- 
erator A is the pair-potential (integral) operator; how 
it acts is specified by the nonlocal kernel A(a;,a;') via 
[kv]{x) = ^ dx' ^{x,x')v{x'). 

Our primary interest is in spectral function p{E) in 
the pseudogap regime as well as the superconducting- 
to-pseudogap transition regime. This quantity can be 
obtained in the usual way from Q{x,x') as follows: 



P{E) 



-TT ^ImTr 



(frg{v,v;u}n)\iu,„^E+iS, (2a) 



r 







dTe*"""g(r,r';T,0), 



(2b) 



where Tr denotes a trace only in the 2x2 particle-hole 
space, (3 = 1/T (i.e., we have chosen units in which 
Boltzmann's constant fee — 1) and the cj„ = (2n -I- 1)-kT 
(with n integral)are fermionic Matsubara frequencies, and 
5 = 0+. We shall assume that the temperature is suffi- 
ciently high to validate the neglect of nonzero Matsubara- 
frequency modes of A, which amounts to treating A as 
a classical (i.e. non-quantal) statistical field. Under this 
static condition, in which the fermion dynamics takes 
place in the presence of an unchanging A field, the spec- 
tral function p{E) may be expressed in the form 



piE)^{piE;A)), 
p{E:A)^y2 5{E-E^), 



(3a) 
(3b) 



where {En} is the collection of energy eigenvalues of 
the following Bogoliubov-de Gennes eigenproblem in the 
presence of a generic classical configuration of A: 



h A 

At ~h 



E 



(4) 



The notation (• • • 
erage over A, i.e. 



denotes the aforementioned static av- 



(5) 



(We do not specifically indicate that this functional av- 
erage is only over static configurations of A.) Thus we 
have expressed the single- particle spectral function p{E) 
in terms of a suitably averaged density of states p{E\ A) 
for a corresponding Bogoliubov-de Gennes eigenproblem. 

Before proceeding with the analysis of the eigenprob- 
lem given in Eq. (0), we address the issue of the form 



of the pairing fluctuations that contribute dominantly 
to the average in Eq. (^). This amounts to a state- 
ment about the physical physical picture of the pseudo- 
gap regime that we are concerned with. Now, as dis- 
cussed shortly before Eq. (^a|), we are considering only 
static configurations of A, and therefore we shall hence- 
forth simplify the notation by writing the pair-potential 
kernel A{x,x') as A{x^x') = A(r, r'). In addition, it is 
convenient to transform A(r, r') to relative and center- 
of-mass coordinates, p and R: 



A(p,R) = A(r,r'), p = r - r' 



2R 



(6) 



It is then convenient to introduce the Fourier transform of 
A with respect to the relative coordinate p, viz., A(k, R): 



A(k,R)^ / d2pe-^'^"A(p,R), 



(7) 



thus obtaining a pair-potential at center-of-mass position 
R and relative momentum k. 

As mentioned in Sec. |[ the scenario for the pseudo- 
gap regime with which we are concerned is based on the 
dominance of configurations of A(k, R) of the form 



A(k,R) = Ao(k,R) e*''^^) 



(8) 



Here, the non-fluctuating factor Ao(k, R) is taken to have 
d-wave form (and can therefore be taken to be real) . 
The fluctuating factor expii?(R) varies slowly with R. 
By assuming this form for A(k, R) we are adopting the 
physical picture of the state as being one in which there 
is no long-range superconducting order, but there are 
local d-wave superconducting correlations, embodied in 
Ao(k,R) [|6|. 

In the following section, we obtain p{E] A) by mak- 
ing use of an elaboration of Andreev's semiclassical ap- 
proach to the Bogoliubov-de Gennes eigenproblem. 
This elaboration is appropriate for the setting at hand, 
viz., one in which there is a strong single-particle poten- 
tial. This scheme was used in Ref. |2^ in order to address 
the low-energy density of states near an extended impu- 
rity in a d-wave superconducting state having a negligibly 
fluctuating, well-formed condensate. In contrast, our fo- 
cus here is on situations in which the fluctuations in the 
amplitude of A are small, but the phase of A is strongly 
fluctuating. Thus, although there is local pairing the sys- 
tem does does not exhibit long-range order. 



III. SEMICLASSICAL APPROACH TO THE 
BOGOLIUBOV-DE GENNES EIGENPROBLEM 



In Sec. H we showed how the single-particle spectral 
function p{E) can be expressed as a density of states 
p{E] A) for the Bogoliubov-de Gennes eigenproblem at 
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arbitrary pair-potential A, averaged over A with a suit- 
able weight, provided A can be treated in the static ap- 
proximation. To make progress with this Bogoliubov- 
de Gennes eigenproblem we invoke, in the present sec- 
tion, a semiclassical approximation under which Eq. (Q) 
reduces to a family of one-dimensional eigenproblems la- 
beled by the classical scattering trajectories of a particle 
in the presence of the single-particle potential V{r). This 
approximation scheme, which was developed in Ref. [ p9[ , 
is valid provided that fcp :$> {Aq, E} for physically rele- 
vant configurations of A (where Ag is the magnitude of 
A), and provided that ^(r) and relevant configurations 
of A vary slowly on the scale of the Fermi wavelength Ap • 




FIG. 1. Schematic illustration of a classical trajectory for 
a quasiparticle scattering from an extended impurity. Also 
shown is a depiction of the background d-wave pair-potential; 
dots indicate incoming and outgoing quasiparticle momenta. 



We now turn to the family of one-dimensional eigen- 
problems arising in our semiclassical scheme. Follow- 
ing Ref. 1^^, it is straightforward to determine that 
the emerging trajectory-dependent eigenproblem has the 
form 



H 



V I \v 

—2ikpds Ao(s) expiz9(s) 

Ao(s) exp — iz9(s) 2ik'pds 

Ao(s) ~ Ao(fcF9sXc(s),Xc(s)). 



(9a) 

(9b) 
(9c) 



Here, i){s) ~ i^(xc(s)), and the parameter s measures the 
position along a particular classical trajectory Xc(s), the 
latter obeying the Newton equation 



fc|a2xe(s) =-VV^(Xc(s)). 



(10) 



Each such classical trajectory (and associated eigenprob- 
lem) is labeled by an impact parameter b and an incoming 
momentum direction n. 



Our next task is to obtain the low-energy eigenvalues 
associated with Eq. ( pa| ) for the case of a d-wave pair- 
potential subject to generic spatial phase variations, i.e., 
Eq. ||. To do this, we perform a local unitary transfor- 
mation of H, i.e., H U H W , where 



U{s) 



-iiJ(s)/2 







g,:^(s)/2 



and thus our eigenproblem acquires the form 



{Ho + Hi 
Ho{s) E 
Hi{s) E 



E 



—2ikpds Ao(s) 
Ao(s) 2ikpd. 

fcp ds^ 
fcp ds'd 



(11) 

(12a) 
(12b) 
(12c) 



The Hamiltonian for this eigenproblem now consists of 
a term arising from the underlying d-wave pair-potential 
(i.e. Hq) as well as a term that contains all the phase- 
variation information (i.e. Hi). Our strategy is to treat 
Hi within perturbation theory, the starting point for 
which is the identification of the eigenstates of Hq. As 
our purpose is to address low-energy states, it is ade- 
quate for us to focus on the zero- energy eigenstate of 
Hq, if any there be. As we shall discuss below, such 
states are guaranteed to arise for certain classical trajec- 
tories Xc(s) 1^,^. The condition for the existence of 
such a state is determined by the asymptotic properties 
of Ao(s): If Ao(s) changes sign an odd number of times 
along an entire trajectory then the Hamiltonian Hq asso- 
ciated with that trajectory has precisely one zero-energy 
eigenstate. Said more formally, if A± = lims^±oo Ao(s) 
then the condition for the existence of the zero-energy 
state is that A+A_ < 0. The explicit form of this eigen- 
state is 



(13a) 



(^±(s) = a± exp±(2fcp)^"'" y ds' Ao{s') 



(13b) 



where ip+ corresponds to the case A+ < (so that 
A_ > 0) and tp_ corresponds to the case A+ > (so 
thatA_ < 0). The pre- factors a± are normalization fac- 
tors, chosen so that (f±{syds = 1. 

Let us emphasize some of the important features of 
these zero-energy states that hold provided the scatterer 
creates only trajectories with at most a single sign change 
in the pair-potential jss). First, the presence or absence 
of these states depends only on the properties of the 
d-wave pair-potential far from any impurity, and is in- 
sensitive to any amplitude variations of the d-wave pair- 
potential that might occur in the vicinity of this impurity. 
Second, the wavefunctions tp±is) exhibit exponential de- 
cay away from the impurity with a decay constant of 
order fcp/|A±|, which is proportional to the BCS corre- 
lation length. As this length-scale is known to be short in 



4 



the cuprates, these states are indeed an extremely local 
probe of local superconducting correlations. 

In the absence of phase fluctuations (i.e. in the pure 
d-wave superconductor), these zero-energy eigenstates 
lead to a sharp peak in the density of states at low- 
energies ll^. To access the impact of order-parameter 
phase variations on this peak, we include the effect of 
Hi on the zero-energy state ip±{s) within quantum- 
mechanical perturbation theory. Thus, the shift in en- 
ergy of the previously zero-energy eigenstate associated 
with a particular trajectory Xc(s) (parametrized by n and 
b) is given by 



ds4(s)i/i(s)V±(s), 



fcp / dsd.s^is)ip±isy 



(14a) 
(14b) 



e(n,6) 



wher e we have used Eq. ( |l3a| ) to get from Eq. (|l4aj) to 
Eq. ( |l4b| ). The spectral function at fixed pair-potential 
p{E, A) consists of contributions from all trajectories 
Xc(s). As we are restricting our attention to low ener- 
gies, it is sufficient to consider only trajectories for which 
A+A_ < (i.e., those for which Hq has a zero-energy 
state). The contribution to the low-energy spectral func- 
tion at constant pair-potential due to the perturbed zero- 
energy states is given by 

p{E, A) ^ fcp 1 ^ JdbS{E- e(n, 6)) 

X (l-sgnA+ sgnA_), (15) 

where the factor (1 — sgnA+sgnA_) ensures that only 
trajectories that satisfy A+A_ < (i.e. exhibit an 
asymptotic sign change in Ao(s)) contribute. 

It is important to emphasize that Eq. (^5|) only in- 
cludes the contribution from the perturbed zero-energy 
states, and that there will be additional, smaller contri- 
butions to the density of states at low energies. Sources 
of these include remnants of the near-nodal quasiparti- 
cles states that exist in pure, d-wave superconductors at 
arbitrarily low energies. 



IV. AVERAGING OVER PHASE 
FLUCTUATIONS 



In Sec. Ill we obtained an expression for the 



Bogoliubov-de Gennes density of states p{E, A) in the 
presence of an extended impurity for the case of a pair- 
potential with local d-wave character and specific real- 
ization of the spatially varying phase i?(r); see Eq. (15). 
In the present section, we calculate the spectral function 
p{E) in the pseudogap regime by averaging p{E, A) over 
suitable configurations of i?(r), i.e., by inserting Eq. ( |l5| ) 
into Eq. (^^. We remind the reader that our scheme 
for computing p{E) applies to settings in which the tem- 
perature is high enough that we may treat the phase 
fiuctuations classically, as discussed after Eq. (pb|). 



A. Gaussian model for phase fluctuations 



In order to evaluate the average over configurations of 
the pair-potential in Eq. ( ^ ) explicitly, we need a model 
for the weight of the various configurations. For the sake 
of simplicity, we choose the phase-field configurations to 
have a (zero mean) Gaussian distribution characterized 
by the correlator (VQi9(r)V/3i?(r')). In terms of this cor- 
relator, it is straightforward to show that the spectral 
function p{E) is a superposition of Gaussian distribu- 
tions, one associated with each classical trajectory on 
which Aq changes sign: 

, fdn /■^^ exp(-£;V2(e(n,6)^)) 
piE) -kp — db ^ ^== (16a) 



2tt 



X (1 - sgnA+ sgnA_) , 
(e(n,6)2)=fc2 Ids Ids' ipl{s)ipl{s'){^{s)^{s')). (16b) 



Here, i9(s) = ds^{s), and the width (e(n,6)2) of each 
Gaussian contribution is sensitive to the degree to which 
local superconducting correlations have been disrupted, 
as can be seen by the pre sence of the phase-phase cor- 
relator (i?(s) i?(s')) in Eq. (16b). Our next task, then, is 



to e valua te the integrals over the parameters s and s' in 



Eq. ( |l6b| ) in order to determine the width (e(n, 5)^) asso- 
ciated with each trajectory (n, b). In the next section, we 
carry out this evaluation within an approximation that is 
valid for the case of long-wavelength phase fluctuations. 



B. Approximate evaluation of trajectory integrals 



In Sec. IV A , we obtained the expression ( |l6aD for the 
spectral function p{E) near an extended impurity in the 
pseudogap regime. In the present section, we make an 
approximation to our expression for p{E) that makes use 
of the local nature of the low-energy states, as well as the 
long-wavelength nature of the pair-potential phase vari- 
ations (appropriate for T ^ Tc). By local we mean that 
the wave functions (p± (s) exhibit exponential decay over 
a length-scale ^ ^ fcF/|A±| (where Aq is the bulk d-wave 
pair-potential), i.e., the BC S co rrelation length; this can 
be seen by examining Eq. ( [l3b| ); see In the cuprate 
superconductors, this length-scale is expected to be much 
shorter than the length-scale for typical pair-potential 
phase variations (i.e. the inter- vortex spacing). Thus one 
has a separation of length-scales: > ^. (Such a sep- 
aration is a natural ingredient of the phase-fluctuation 
picture of the pseudogap regime because for inter- vortex 
spacings on the order of ^ the conventional meaning of 
local Cooper pairs to break down.) 

To make the approximation to our expression for p{E), 
consider the integrations over the trajectory parameters 
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s and s' in Eq. (16b). Now, the correlation function 
(i?(s) "dis')) varies appreciably only over length-scales on 
the order of ^,9 or longer, whereas the wavcfunctions ip± 
decay exponentially, as mentioned, on the length-scale 
^. Hence, one can make an asymptotic approximation 
to the s and s' integrations which amounts to pulling 
the correlator out of the integrals. Thus, owing to the 
normalization of ip±, we have 



(17) 



independent of n and b. 

Next, we turn to the interpretation of the quantity 
(i9(0) 1^(0)). The derivative of the phase along a particu- 
lar trajectory is given by 



l?(s)=Xe(s)-Vl?(x)| 



(18) 



By inserting Eq. ( |18| ) into Eq. (17), we see that the cor- 
relator of interest is (Va'?(r)V/3i?(r)) which, by spatial 
isotropy, has the form (5„/3 (| Vi9(r)p)/2. According to 
the phase-fluctuation scenario, vortex excitations are the 
dominant mechanism for generating phase gradients . 

How does this information about the phase correlations 
translate into information about the spectral function? 
By noting that trajectories Xc(s) involving sign-changes 
in the d-wave pair-potential obey |xc '--^ 1 and using the 
aforementioned phase-gradient correlator, we arrive at 
the width 



4(|V^(r)p)/2 



(19) 



and, hence, our final expression for the spectral function 
p{E): 



P{E) 



exp(-£;V2(e2)) 



(20) 



xfci 



'/£ I '^^ (l-sgnA+sgnA_). 



The first factor is a Gaussian in the energy E with a 
linewidth (e^)^/^ associated with the r.m.s. fluctuations 
in the phase gradient. The second factor (i.e. the in- 
tegral over the impact parameter h and the incoming 
momentum direction n) determines the scale for p{E), 
essentially by counting the number of sign-changing tra- 
jectories, and is expected to be only weakly temperature 
dependent. In the next section, we compute the linewidth 
within a Berezinskii-Kosterlitz-Thouless-like model of 
phase fluctuations near Tc. 



C. Linewidth near the superconducting phase 
boundary 

To make further progress, we now attempt to calculate 
the spectral function linewidth due to the phase fluc- 
tuations accompanying the destruction of superconduct- 
ing order in the neighborhood of T^. We shall do this 



by choosing a particular weight for the phase fluctua- 
tions, viz., that associated with the two-dimensional XY- 
model ]40|-|4^. Our analysis is reminiscent of that due 
to Franz and Millis ||ll[ who addressed the bulk single- 
particle spectral function in the pseudogap regime. In 
the present context, this analysis is provided solely for 
illustrative purposes, and is only meant only to provide 
a rough estimate of the linewidth. 

Let us consider the XY-vcioAel action. 



K 
1 



(21) 



where K = ps (T) /T in which Ps (T) is the temperature- 
dependent superfluid density ||43|| . Note that we are 
not suggesting that the true critical fluctuations of the 
cuprate superconductors necessarily lie in the universal- 
ity class of the two-dimensional XK-model, but simply 
that the intermediate length-scale fluctuations proposed 
as leading to pseudogap phenomena may adequately 
be modeled by Eq. (|2l|). The form for (|Vi?(r)p) for 
the XF-model may be calculated following the Debye- 
Hiickel-type analysis of Halperin and Nelson 1 45 1 , giving 



|Vi9(r) 



(22) 



in which K* is the short length-scale stiffness (obtained 
using the Kosterlitz renormalization-group equations), A 
is a short-distance cutoff, and is a length-scale char- 
acterizing the typical inter- vortex spacing. The principal 
temperature dependence in Eq. ( p^ ) arises via near 
Tbkt [i-e. the transition temperature of the model (|2l|), 
which is expected to lie not far below Tc] is propor- 
tional to exp ^yQ/{T — Tbkt), where is a constant of 
order unity (which we do not try to calculate in detail). 




FIG. 2. Spectral linewidth versus (T 
Immediate vicinity of the transition. 



rBKT)/e. Inset: 



The proliferation of unbound vortex excitations upon 
warming through the BKT transition is reflected by this 
strong diminution of and causes a concomitant dra- 
matic increase in the linewidth of the spectral function: 



(e2)i/2 ocexp 



BKT, 



(23) 
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Observation of a fluctuation-broadened peak in the spec- 
tral function, e.g., via STS measurements, would pro- 
vide striking evidence in support of the phase-fluctuation 
scenario. Moreover, the temperature dependence of the 
linewidth would, e.g., provide access to the details of the 
vortex-unbinding transition. 



V. NUMERICAL ESTIMATE OF THE 
LINEWIDTH; CONCLUDING REMARKS 

In the previous section we saw how one could esti- 
mate the temperature dependence of the linewidth of 
the spectral function near the superconducting transi- 
tion temperature Tc- In the present section we make a 
rough numerical estimate of this linewidth at one partic- 
ular temperature in the pseudogap regime by appealing 
to the data obtained in the high-frequency conductivity 
experiments of Corson et al. Q on Bi2Sr2CaCu208+5. 
We shall be speciflcally concerned with the interpreta- 
tion of these data inasmuch as they provide access to the 
characteristic vortex density. 

Until now we have been working with a system of units 
in which fi^ /2m — 1. Restoring conventional units in 
Eq. (n9) gives for the linewidth 



2N1/2 



2m 



fcF(|V,9(r)|V2)i/'. 



(24) 



To estimate this width, we turn to the Corson et al. ex- 
periments and its analysis by Corson et al., which is 
based on the notion that, at sufficiently high frequencies, 
the conductivity probes short-length-scale pairing corre- 
lations, and leads to an estimate for the characteristic 
density of free vortices rtf. Assuming that is is vortex 
excitations that lead to phase fluctuations, one expects 
that, up to a constant of order unity, (| Vi9(r)p) ~ nf. 

Now, Corson et al. obtain values of n/ that are on 
the order of 0.003 Ovc (for T ^ 75K), where Cvc is the 
area of the core of a vortex. If we take the vortex core 
to be a disk of radius — 1 nm, this leads to the value 
rif ~ 10'^/J,m~2. Then, using the order-of-magnitude es- 
timate fcp ~ lnm~\ Eq. (|^ gives {e^^^^ ~ 9meV. 

The value of this estimate is that it shows that, for 
at least one cuprate material, there is a temperature 
at which the free-vortex density is small enough that 
phase fluctuations only weakly perturb the energies of 
the quasiparticle states. Hence, the linewidth arising 
from perturbed zero-energy states can be rather smaller 
than the scale of the superconducting energy gap (and 
hence small enough to justify the our picture of per- 
turbed zero-modes) but be large enough to be resolv- 
able in STS measurements, such as those of Yazdani et 
al. j27| . Larger densities of free vortices, and hence large 
line- widths, would result from working at higher tem- 
peratures. Moreover, near Tc the vortex density is ex- 
pected to show a very strong dependence on tempera- 
ture, which should confer a strong temperature depen- 
dence on the linewidth. Of course, if conducted at the 



temperature of the Corson et al. experiments, thermal 
broadening would complicate the task of accessing the 
intrinsic linewidth (i.e. the linewidth due to phase fluc- 
tuations); recall that 10 K is equivalent to 1 meV. Thus, 
as emphasized in Sec. ^, one should consider performing 
experiments on materials having a lower Tc, so that the 
pseudogap regime can be explored at temperatures at 
which thermal broadening is less significant. These con- 
siderations indicate that it is at least conceivable that 
STS experiments near extended scatterers could provide 
a sharp test of the phase fluctuation scenario for the pseu- 
dogap regime. 
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